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Abstract: We study chiral symmetry restoration with increasing temperature and den- 
sity in gauge theories subject to mutuahy perpendicular electric and magnetic fields using 



cn i holography. We determine the chiral symmetry breaking phase structure of the D3/D7 and 



D3/D5 systems in the temperature-density-electric field directions. A magnetic field may 



lO ' break the chiral symmetry and an additional electric field induces Ohm and Hall currents 

C^^ ' as well as restoring the chiral symmetry. At zero temperature the D3/D5 system displays 

^-^ ■ a line of holographic BKT phase transitions in the density-electric field plane, while the 

D3/D7 system shows a mean-field phase transition. At intermediate temperatures, the 

transitions in the density-electric field plane are of first order at low density, transforming 

K> , to second order at critical points as density rises. At high temperature the transition is 

H I only ever first order. 
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1 Introduction 

There has been considerable interest recently in the chiral symmetry breaking/restoration 
phase structure of gauge theories with quarks [1] that can be described holographically 
[2-4]. The simplest model to study is the 3+1 dimensional M = 1 gauge theory made from 
adding quark hypermultiplets to A/" = 4 super Yang-Mills theory [5, 6]. This theory can 
be described holographically, in the quenched limit, by probe D7-branes in an AdS5 x S° 
geometry [7]. Chiral symmetry breaking can be induced by the application of a magnetic 
field which breaks the conformal invariance [8]. Temperature, T, is introduced by an 
AdS-Schwarzschild black hole geometry [4] and density, d, through a non-zero temporal 
component of a gauge field on the Z?7-brane world-volume dual to a C/(l) baryon number 
background gauge field [9-12]. The chiral symmetry restoration transition is first order 
[13, 14] with temperature and second order with density [15, 16]. There is a critical point 
linking these behaviours in the T-d plane [15]. A recent study [17] has also looked at the 
chiral symmetry breaking phase structure in a background geometry with a running gauge 
coupling displaying first order transitions as well as a baryonic phase. These results are 
potentially of use in describing 3+1-dimensional gauge theories that break chiral symmetry 
including, of course, QCD. 

The D3/D5 system with a two dimensional defect is also of interest since it describes 
a 2+1-dimensional theory of quark hypermultiplets in the background of 3+1-dimensional 



AA = 4 adjoint matter [18-23]. This may be relevant to condensed matter systems. Again 
a background magnetic field induces chiral symmetry breaking [22, 23]. Here the novel 
discovery at zero temperature but finite density was a holographic BKT transition across 
which the quark condensate grows exponentially with the distance from the critical density. 
This behaviour results from the fact that density and magnetic field are both dimension 2 in 
2+1-dimension and so compete in the deep IR as explained in [22]. Any finite temperature 
converts the transition to a second order mean field transition. Examples of cases in which 
a mean field transition interpolates to the BKT point through non-mean field second order 
transitions was also explored in [24]. 

The rich structure of these results is very interesting and it is natural to ask how 
generic the phase structures found are. In part to address this question in [25] the authors 
introduced an electric field parallel to the magnetic field already studied. The electric 
field simply adds another parameter acting against the chiral phase (the E field tries to 
dissociate the mesons). It is of course also interesting in a condensed matter context to 
study the effects of E fields. The chiral phase structure in the T-E plane was very similar 
to that in the T-d plane. The E-d structure of the theory was rather singular though with 
density not being generated no matter how large a chemical potential was applied. The 
DEI analysis was suspect at the origin of the space where the ^(1) baryon number gauge 
field diverged ^. 

In this paper we wish to add to the investigation of the effects of an E field by study- 
ing the case of an electric field perpendicular to the magnetic field. This case does not 
show any signs of singularities in the DEI analysis but is also technically more difficult 
because of the presence of a Hall current. It is however possible to extract all the currents 
analytically [26, 27]. We will work in the Canonical Ensemble where we treat density as a 
free parameter. The crucial transitions on the gravity side are then associated with D7 or 
Z?5-brane embeddings. If the branes lie flat they describe a chirally symmetric phase. If 
they bend in the space transverse to the D3 branes they describe chiral symmetry breaking 
configurations [28]. 

Our results are summarized in Fig 7 where we show the E-T-d phase space at fixed 
magnetic field. In the D3/D7 case the chiral transitions are first order at low density whilst 
second order at high density. We have identified a critical line between these transition 
types. The D3/D5 system is similar except that it has in addition a line of EKT transitions 
in the E-d plane at zero temperature. The phase boundaries on the E-d surface form a 
square - we show that the EKT transition is driven by deep IR effects and that they 
are independent of E until the E field itself is super-critical when it overcomes the EKT 
transition entirely, explaining this shape. In this case there is no interpolation from mean 
field to a EKT transition. 



^In the case of a parallel E and B field there is a Wess-Zumino term in the D7 action which is proportional 
to f d'^xE ■ B. This is a total derivative and hence if the background E and B fields fall of sufficiently 
quickly asymptotically this gives the zero winding number of the (7(1) configuration. Presumably these 
terms can be neglected as we did in [25]. Were the fields not to fall off as quickly as 1/r^ there might be 
additional terms that would effect the analysis. In the case of perpendicular E and B fields we study here 
there are no such Wess-Zumino terms. 



2 D-brane setup 

We study a single probe flavour brane embedded in the Poincare patch of the Schwarzschild 
AdS^ X 5^ geometry, generated by a single stack of A'^c S> 1 black D3-branes. We choose 
the following metric parametrization: 



ds' = GMNdx'Ux'' = ^ i-fdt' + dx') + ^^ + R'dnl , f{r) = 1 - ^ , (2.1) 



J.4 



where R is the radius of the AdS^ and S^ spaces, dx^ = dx"^ -\- dy^ + dz^ and df^s is the 
unit radius S^ . r goes from the horizon dX r = r^ to the boundary at oo. The temperature 
of the black hole and the associated boundary field theory are both given by: 

r„ = ^. (2.2) 

We will specialize to the cases of Dg-probe brane phenomenology with q = 5 and 7. In the 
q = 5 case the D3 and D5 have 2+1 coincident dimensions and the description corresponds 
to a theory with fundamental matter living on a 2 + 1-dimensional defect of the AD CFT 
and may be used to study the phenomenology of some condensed matter systems [18-23]. 
q = 7, when the D3 and D7 have 3+1 coincident dimensions, is the canonical model for a 
flavoured 3+ 1-dimensional theory and, with the inclusion of a conformality breaking scale 
and can be considered as a QCD-like model [7, 8, 13-16]. 

The Dg-brane action is given in terms of the metric pulled back onto its worldvolume. 
This includes the scalar fluctuations of its transverse dimensions and to this we add the 
gauge field contribution. All fermionic fields will be turned off in what follows. We will 
also not include the Wess-Zumino term which vanishes for all cases of interest. The DBI 
action is thus: 

Sbbi = -ToqNf J d^+^C e- V-det(7mn) , (2.3) 

where Tpg = (27r)~''a'~^'^~^^)'^ and 7 denote the pullback metric plus gauge field 

^mn = d^X^'dnX^GMN + 2T,a' F^n ■ (2.4) 

We will consider different choices for the gauge field on the Dg-brane, but the one-form 
gauge potential {A^ such that in equation (2.4) F = dA) which captures all phenomena 
discussed here is [8, 10, 14, 26]^ 

2Tra'A = At{p)dt + {-E^t + A^{p))dx + {B^x + Ay{p)) dy , 

where the Ai^ E^, and B^ variables include the factor of lira' and correspond to the 
background field. At introduces finite baryon density, E^ is a constant background electric 
field and B^ a magnetic field in the x and z directions respectively. We will drop the x and 
z indices from now on. A^ and Ay will encode the Ohm and Hall currents generated by 



^In [29] an isospin chemical potential was turned on using a pair of D7-branes with a U{2) gauge field 
and the phase structure of this system studied in detail. 



the electric field and magnetic field. This choice of the gauge potential is sufficient to find 
the vacuum Dl and Df> configurations with our E and B field choices (though if we were 
interested in mesonic states described by fluctuations of the gauge field degrees of freedom 
we would need a more complete ansatz). In the case of the D5-brane, the magnetic field is 
a scalar, while in the D7-brane case it is a pseudovector, pointing in a direction transverse 
to the electric field. The fact that the electric and magnetic fields are pointing in different 
directions results in a vanishing Wess-Zumino term. 

To describe the embeddings in a convenient way we parameterize the S^ as 

d^l = dO"^ + cos^ em1^__-^y^ + sin^ 0<^^%-q)/2 > (2-5) 

such that the probe brane wraps the ^ — sphere and its transverse directions are parametrized 
by the -^-sphere and the 9 direction. We will find it useful to make the coordinate trans- 
formation 



dr^ dnP' 



^ = -^ -^ V2w = ^r^ + ^r^-rj,, (2.6) 

with V2wh = fH- This change makes the presence of a flat 6-plane perpendicular to the 
horizon manifest. We will then write the coordinates in that plane as p and L according to 



w = \/ p^ + L"^ ^ p = w cos 9 , L = w sin 9 , (2.7) 

The metric is then 

d^' = ^(-Sttdt' + fcdf2) + ^(dp2 + p^^nl + dL2 + L'dnj) , (2.8) 



where 



/ 4 4 \2 4,4 

[W — Wjj) _ W + Wjj 



- _ \^ ^H) - _ "^ ^ '^H ^9 QN 

y** ~ ...4/, ..4 , ...4 N ' y^JX — _,_4 • \'^-^} 



In this coordinate system, we study the embedding scalar L = L{p) assuming all other 
scalars to be constant (and indeed zero by symmetry). Therefore, the DBI action at zeroth 
order in fluctuations (2.3) effectively describes the dynamics of a single scalar fleld plus a 
gauge field as a function of a single radial direction 

Sbbi=M fdtdxfdpC[L{p),L'{p),A'^;p] , (2.10) 

where TV = NfToqQ'^ J dO(g_i)/2. 



2.1 The Currents 

We will begin by determining the currents induced by the background E and B fields. 
Since A^i occurs only through its derivative we define the associated conserved current as 

^" = 1^ ■ <"^) 

and from this we can define the Legendre-transformed Lagrangian 

sgn(g) 



C^T^C- J^^A'^ = - ^' VC^[^a^jGpp + GllL'^ , (2.12) 



e = 


~GttGxx ~ GttB — GxxE , 


x = 


9-1 9-1 

—d Gtt — GttGxx G^l^ — Gxx{Jx + "Jy) 


a = 


—dGttB + GxxJyB . 



where d = Jt and 

(2.13) 

(2.14) 
(2.15) 

In order to fix the values of the currents for a given external field configuration we 
have to distinguish between the finite and zero temperature configurations. 

2.1.1 Finite temperature 

With finite electric field and at finite temperature there is always a position, Wg > 0, where 
^ changes its sign, which we call the singular shell, [14]. In order to keep the Lagrangian 
real we have to require that both x and a change their sign at r^. i.e. x(^s) = cl{ws) = 0. 
The logic is that a must vanish at the point where ^ vanishes to keep the action real since 
a^ enters with a negative sign. Further, so that the action is real in the regime where ^ is 
negative x must also have switched sign so that ^x is positive. This constraint of reality is 
enough to determine the currents completely from the background field values [26] 



'- = - e + T'Jm "• <^-'«' 

J. = - ..lf..r„ E. (2.17) 

where T = ttRTh, Ps is the position at which the embedding touches the singular shell, 
and 



I'^rPfa 1 r^) 1 4(g+T4)(g+T4+2B2)p(P-i) 


V -"^ ^^ ' ^ ' GWG'^T^ 


g + T^ + 2^2 

2dB 

Q + T^ + 2B^ ' 



g = -B^ + E^ + ^J-AB'^E'^ + (T^ + 52 + ^2)2 , (2.18) 

We see therefore that when we have finite magnetic field, electric field and baryon 
density at finite temperature there is both an Ohm current, J^, and a Hall current, Jy, 
caused by the Ohm current moving transverse to the magnetic field. The Ohm current has 
two contributions - one is proportional to d and represents the quark density in the vacuum 
reacting to the electric field. The second part of the current is due to the pair-creation of 



quarks and anti-quarks from the vacuum. Note that the Hall current only has a non-zero 
contribution from the quark density. This just reflects the fact that the Hall currents for 
the pair created positive and negative charges cancel out completely. 

2.1.2 Zero temperature 

Unlike at finite temperature [26, 27], at zero temperature there are two cases which have 
to be studied separately: 

E > B: When the electric field is larger than the magnetic field there is a singular 
shell at finite radius, defined via ^ = 0, just as in the finite temperature case. So, from 
X = a = 0, 



*Jx 



\ 




, -^2/ = -^' (2-19) 



where we have used pi = ^/ E"^ — B'^, since the embedding is flat in this case. 

E < B : When the electric held is smaller than the magnetic held, the equation ^ = 
does not yield a flnite radius singular shell (r^ = 0), so the reality condition is not ^ = a = 
at /3 = 0, but S,x ^ ^^ foi' 8-11 P- 

iX-a"^ ((p2 + L^f + B^- E^){p^ + dF- Jl - J2) - {dB + JyEf (2.20) 

>-Jl{B'^-E^)-{JyB + dEf. (2.21) 

Thus the condition ^x ~ ^^ ^ implies that 



dE 



Jx = 0, Jy = -—. (2.22) 



Note that .^ = a = Oatps = does not yield a consistent result. This result holds for both 
the Dl and D5-branes. That J^ here is zero but the Hall current is non-zero is a standard 
steady state Hall effect result at large magnetic field. J^ = because the Lorentz force due 
to B and the electric force due to E at large B are equal and opposite. However, if E > B, 
the force due to E starts dominating over the Lorentz force and a finite Ohm current is 
generated (2.19). 

If we consider J^ as an "order parameter" of the quantum phase transition (nonzero at 
E > B and zero otherwise) then there is a second order phase transition at E = B whose 
critical exponents are, from (2.19), 

4 ~ (£; - 5)1/2 ^ p5) ^ (2.23) 

4~(E-S)i/2 {D7,d^0), (2.24) 

J^r^{E-Bf/^ {D7,d = 0). (2.25) 

Without density, the critical exponents are understood on dimensional ground (1/2 for D5 
and 3/4 for D7). At finite density, the phase transition is governed by a density effect 
(1/2 for both D5 and D7). For the Z?5-brane case, they happen to be the same. Note 



that a thermal insulator-conductor transition is of first order as shown in section 3.3'^ In 
the following sections, we will show that these critical exponents remain also valid for the 
chiral condensate. 

2.2 The embedding and chiral condensate 

Having found the currents it now only remains for us to find the D7 or L>5-brane embedding 
in the geometry. The U{\) chiral (R-) symmetry of the gauge theory is dual to the U{\) 
rotational symmetry in the ^(9_q)/2 of (2.5). This rotational symmetry is broken when 
the embedding is not L = 0. The equation of motion of L{p) is obtained by varying (2.12) 
after plugging in (2.16) and (2.17) (or (2.19), (2.22)). In general the classical embedding 
behaves at large p (UV) with two parameters m, c 

Lp^^ ~ m + j^ . (2.26) 

We are most interested in studying the phenomenon of spontaneous chiral symmetry break- 
ing and thus require the symmetry to be restored in the UV. For this we must impose that 
m = and determine the value of c by demanding IR regularity of the embedding solution. 
If the regular (and lowest free energy solution) gives c = then the system is in a chirally 
symmetric phase, whereas for c 7^ the chiral symmetry has been spontaneously broken. 
At zero temperature and density we impose the IR boundary condition L'(0) = shooting 
out to find regular solutions with tti = 0. At zero temperature but non-zero density we seek 
solutions that end in the IR at the origin so -L(O) = 0. Finally at non-zero temperature we 
shoot out from points on the singular shell and seek solutions that tend to L = at large 
p. Points where there are phase transitions can be identified by plotting the parameter c 
as the background fields are varied. For more technical details we refer to [15, 23, 25]. 

3 The D3/D5 and D3/D7 phase diagrams 

Here we study the phase diagram of the above D5 and D7-brane systems in the d, E, B, T 
directions. The two brane systems will be discussed in parallel, even though we will see that 
there are some stark differences between the phase structures. We will always be studying 
the cases of massless fundamental matter and thus there are four dimensionful parameters. 
Because the underlying theory is conformal, we can use one of these parameters to set 
the scale for all other parameters. We chose B = 1 and set all other scales by this value. 
Rather than giving the full three-dimensional phase diagram, to start with we will focus 
on two dimensional sections of it, contrasting the different behaviours in the two brane 
systems. 



^Note that the zero-temperature insulator conductor transition can also be first order if we include finite 
mass for the fundamental matter. In [30] it was shown that at zero temperature and zero magnetic field, 
the transition from Minkowski, to singular shell embeddings has a region of first order behaviour as well as 
second order behaviour. 
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Figure 1. T-d phase diagrams at _B — 1 for the D5 (a) and D7 (b) probes. The red hncs correspond 
to second order phase transitions and the blue to first order phase transitions. There is an isolated 
BKT phase transition in the D5-hra.ne case at exactly T = due to the competition between B 
and d which in 2 + 1 dimensions have the same scaling dimension. The BKT phase transition is 
absent in the D7 case because the dimensions are [d] = 3 and [B] ~ 2 which means that such a 
competition cannot arise and thus the IR phenomenology is dominated by d. 



3.1 T-d phase diagram at finite magnetic field 

The phase diagram at finite chemical potential, temperature and magnetic field has been 
studied in detail in [15, 16, 22, 23]. In the following we focus on the canonical ensemble, 
the phase diagram of which is shown in Fig 1. 

The finite magnetic field increases the binding energy of the fundamental charges and 
thus stabilizes the mesons, whereas the thermal fluctuations at finite temperature act to 
destabilize them. The finite charge density means that one can never have fully stable 
bound states and thus the width of all mesonic fluctuations is finite. The transition is 
thus not a meson melting phase transition (as is the case at zero baryon density) but a 
chiral symmetry breaking phase transition. For large baryon densities or temperatures, 
the embedding is always the trivial embedding but the transition from the trivial to the 
non-trivial embedding may be first order or second order in the L)7-brane case, with the 
addition of a BKT phase transition at T = in the D5 brane case^. This has already been 
discussed extensively in [22-24]^. 

3.2 The phase diagram with finite E and B 

The inclusion of both finite electric and magnetic field has been discussed previously in 
[25] however in that case the electric and magnetic fields were parallel and thus there was 
no Hall current. In the current paper we study the case where the two are perpendicular, 
giving a very different phenomenology and a somewhat different phase structure. In the 
following sections we will deal with various slices in the finite E part of the phase diagram. 



*At high temperature the first order phase transition line turns out to be unstable thermodynamically, 
which can be addressed in the grand canonical ensemble more properly [15]. However, here we deal only 
with the canonical ensemble. 

''A BKT transition in the back-reacted U{2) Einstein- Yang-Mills theory in 4-fl dimensions has been 
discussed in [29]. 




(a) Embeddings with m — 0. Col- (b) c vs E. Red dots are numerical 
ors match large dots of (b). values and the back curve is c = 

0.58(1 - £;^)i/2 
Figure 2. Second order phase transition for the D3/D5 system with B = 1 varying EatT = d = 0. 





(a) Embeddings with m — 0. Col- (b) c vs E. Red dots are numerical 
ors match large dots of (b). values and the back curve is c = 

0.225(1 - £^)=''''* 

Figure 3. Second order phase transition for the D3/D7 system with B — I varying ii^ at T = rf = 0. 



The phenomenology of branes with finite electric fields has been covered extensively 
in [14, 31, 32] and in [27] the interpretation of an open-string horizon with a membrane 
interpretation was elucidated. 

3.2.1 Transition on the E-axis: T = d = 

In section 3.1 we discussed the T-d phase diagram. By looking at the behaviour on the E- 
axis (again, at finite, fixed B), that is at T = d = we can get an idea of the full 3D phase 
structure. Taking the T = d = case of the Legendre transformed Lagrangian,(2.12), and 
specializing to the E'-axis we get the simplified Lagrangian: 



LLT = p(^-^)/^ 



1 + 



1-^2 

(T+W 



x/OTl^ 



(3.1) 



where we are restricting to the case where E < 1 (i.e. E < B) such that there is no singular 
shell and no current induced (2.22). For E > 1 the embedding is always flat and thus there 
is no condensate whereas for E < 1 there is a finite condensate which vanishes as S — )• 1. 
In this case the critical exponent of the system can be computed analytically. If we 
substitute 1 — E^ for B, then the actions are identical to the case at finite B and T = E = 
formally. This is nothing but a Lorentz transformation. On dimensional grounds the 
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Figure 4. T-E phase diagram at B = 1 and d = Q. All fist order transitions(blue) except one 
second order point (red) at T = 



condensate must be related to the electric field as: 



(3.2) 



The proportionality constant is determined by numerical analysis at any one point of 
They turn out to be 



c = 0.580(1 - -^2)^2 (^5) ^ 
c = 0.225(1 - E^f/^ {D7) , 



(3.3) 
(3.4) 



where Ec = 1- The D3/D5 system has the mean-field exponent 1/2, while that of the 
D3/D7 shows the non-mean-field exponent 3/4. This is in agreement with the Ohm cur- 
rents (2.23) and (2.25) which are also understood on dimensional grounds. 

We are able to check our numerical procedures against the analytic results and this is 
shown in Fig 2 and Fig 3. In each of these plots, (a) shows the embedding configuration 
for a number of different values of the electric field strength, E, while (6) gives a plot of 
the condensate as a function of E going from to 1. The large coloured dots in the (6) 
plots correspond to the three embedding configurations in (a). The black curves ((3.3) and 
(3.4)) agree with the red dots, which are numerical values obtained for various E. 

3.2.2 T-E plane 

The Lagrangian (2.12) with d = is rather complicated and not terribly illuminating so 
we don't present it in full here. Both electric field and finite temperature act to destabilize 
the fiuctuations of the fundamental bilinear operator. The critical line in the T — E plane, 
shown in Fig 4, corresponds to a first order phase transition between a chirally symmetric 
phase and a chiral symmetry broken phase. This first order transition becomes a second 
order transition at T=0 which is of mean field type in the D5 case and a non-mean field 
transition (with critical exponent 3/4) in the D7 case. 
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Figure 5. The d-E phase diagram for the D5 (left) and D7 (right) systems. Chiral symmetry is 
broken inside the enclosed region. B = 1 and T = 0. 



3.2.3 d-E plane 

At zero temperature the Lagrangian (2.12) simplifies greatly, with the currents given in 
(2.19). 



^LT 



^Mi-|^)+p''-Mi + 



B'^-E^ 



] Vl + L'^ ■ 



(3.5) 



(p2+L2)2 

Although the Lagrangian is simplified, the numerics at T = are somewhat more compli- 
cated than the finite temperature case. At T = we are not able to integrate the equations 
of motion from the IR at r = and instead have to solve from the UV. At finite tempera- 
ture the singular nature of the equation of motion at the horizon means that given a single 
boundary condition plus the condition of regularity we can shoot out to the UV and pick 
off the solution which has zero mass. Shooting from the UV into the IR however we choose 
the zero mass solution and then have to find the value of the condensate such that the 
brane embedding hits the point r = 0. Because of the singular nature of the equation of 
motion at this point, the numerics are somewhat sensitive to the accuracy which we re- 
quire. The second subtlety is that because there are actually an infinity of solutions which 
hit the singularity (associated with the self-similar spiral structure in the c — m plane, [8]) 
we must be careful that we pick the one associated with the lowest free energy. 

Having tamed the numerics we are able to study the phase transition from the chirally 
symmetric phase to the chiral symmetry breaking phase and this is plotted in Fig 5. For 
the D5 case the analysis of the condensate as a function of d and E is shown in Fig 6. The 
phase transition as d is increased is independent of the size of the electric field (for E < 1) 
and is of BKT type. As E increases for d < dc = \/7 the phase transition is of mean field 
type and again its position is independent of the value of d. This should be contrasted 
with the D7 system where the second order phase transition takes a non-trivial profile in 
the d — E plane. The D7-brane analysis is relatively simple and so is not presented here. 

The line of BKT transitions in the D3/D5 system and its -E-independence can be 
understood from the IR dynamics. After scaling all quantities by B we have the embedding 
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(a) c{E,d). BKT transitions on the line 
d = \/7 for (0 < £ < 1). Mean field second 
order transition on the line E = 1 for (0 < 
d< ^/7). 
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(c) c YS d cross-section near phase transition d = 
\/7 for a fixed E. E = 0.1, 0.2, • • • ,0.9 from or- 
ange to purple. BKT scalings are shown, the 
BKT transition point dc = \/7 is independent of 
E as shown with the IR analysis. 




(b) c YS d cross section for a fixed d. The 
critical exponent is 0.5 ±0.0001 near phase 
transition point E = 1 . 
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(d) Numerical confirmation of BKT scal- 

-■n I ^ + '^\ 

ing. c ~ e V ''s-'^ . The slopes are 
1± 0.001. 



Figure 6. The analysis of the BKT and mean field phase transition for the D3/D5 system. B — 1 
and T = 0. 



Lagrangian, at small p, 

Clt ~ \/l - E'^^/d? + p\p^ + L2) Vl + L'^ 



(3.6) 



Note that the E part is factored out from the Legendre transformed action so that the 
dynamics is identical to the E = Q case. The slipping mode for a small fluctuation of the 
embedding L in the IR p ^ (l;*^) behaves as a scalar in AdS2 space with mass — ^^^^ . 
The Breitenlohner-Freedman (BF) bound for a scalar in AdS2 is —1/4, so below d = \fl 
the BF bound is violated and the phase transition to a broken symmetry phase occurs. 
Note that we get a good fit for the whole of the region close to the phase transition with 



c ~ 2.6e V7-<*^ \/l — E. The exponent comes from the analytic calculation given in [22]. 

3.3 T-d-E phase diagram 

Finally, in Fig 7 we combine the above information with the phase structure outside the 
two dimensional slices we have analyzed above. The critical point linking the first and 



12 




0.3 

(a) D3/D5 (b) D3/D7 

Figure 7. 3D phase diagram of D3/D7 and D3/D5 system at S = 1 as a function of T,d, and E. 

second order phase transitions give a line on the surface of the phase transition face. This 
links surfaces of second order and first order phase transitions. The line of BKT phase 
transitions is an isolated region in the D3/D5 diagram and the rest of the diagram is 
either of mean-field second order type or of first order type. Fig 7 summarizes our full 
analysis. 

4 Summary 

We have studied the D3/D7 and D3/D5 systems in the presence of mutually perpendicular 
electric and magnetic fields to deduce the temperature and density phase structure of the 
dual gauge theories. Those theories are the 3+1 dimensional A/" = 4 super Yang Mills 
theory with quark hypermultiplets introduced in all field theory directions and on a 2+1 
dimensional defect respectively. The full results are displayed, for constant magnetic field, 
in Fig 7. There are two phases. At large d, T, E the theory is in a chirally symmetric phase, 
while at low d, T, E the phase displays chiral symmetry breaking through a non-zero quark 
condensate. In the gravity description these transitions are associated with embeddings 
of the probe D7 or D5-brane that either lie flat (the symmetric phase) or are given by 
non-trivial profiles in the holographic directions breaking the symmetry. 

At zero temperature the phase transition is continuous for both D3/D7 and D3/D5 
systems. However, the D3/D5 system has a special feature that there is a holographic 
BKT transition for low E field values. In this regime the deep IR is dominated by both 
the magnetic field and density regardless of E. When the electric field becomes sufficiently 
large to trigger the transition on its own the transition becomes second order and mean 
field. At intermediate temperature, at low density the transition is first order whilst at 
high density it is mean-field second order. The critical lines between these regimes can be 
inferred from the Fig 7. At high temperature the transition becomes all first order. 
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If we were to include finite mass, the structure may fiave even more subtleties within. 
The more features we can find in our phase diagrams, the more chance there is we can link 
such basic holographic models with real world systems. We leave the extensions of this 
phase space for future work. 
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